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Abstract 

J> \ The operator product expansion (OPE) for the difference of vector and axial current 

correlators is analyzed for complex values of momentum q 2 . The vector and axial 
. spectral functions, taken from hadronic r-decay data, are treated with the help of Borel, 

O ! Gaussian and spectral moments sum rules. The range of applicability, advantages and 

disadvantages of each type are discussed. The general features of OPE are confirmed 
■ by the data. The vacuum expectation values of dimension 6 and 8 operators are found 

^ ! to be O e = -(6.8 ± 2.1) x 1(T 3 GeV 6 , O s = (7 ± 4) x 1(T 3 GeV 8 . 

a: 

§. 1 Introduction 

& : 

Precise measurements of vector V and axial A spectral functions in r-decay have been 
recently performed by ALEPH |T| and OPAL |2| collaborations. Define the polarization 
^ | operators of hadronic currents: 

KM = i j e%qX < TU,{x)U v {^ >dx= (q,q u - g, v q 2 ) IlgV) + ^4V) , W 
where U = V, A ; V ll = u^^d , A^ = uj^d . 



The imaginary parts of the correlators are the so-called spectral functions (s = q 2 

W/A 



vx/axis) = 27rlmn^ (s + iO) , a (s) = 2nlmUf (s + iO) . (2) 



which have been measured from hadronic r-decays for < s < m 2 . The spin-0 axial spectral 
function a (s) is basically saturated by r — > ixv T channel, which gives 5-function. It will not 
be considered here. 

In this paper the experimental data for v\ — a\ will be used to determine numerical values 
of the quark condensates in QCD. An early attempt to realize such programm was performed 
by Eidelman, Vainstein and Kurdadze || using e + e~ annihilation data, but the experimental 
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errors were rather large and the result not very conclusive. Also, higher order condensates 
and higher order perturbative corrections were not included in the analysis. More recent 
analysis M of e + e~ annihilation data demonstrates, that the spread of the values of the 
quark and gluon condensates is larger than found in || . Therefore the consideration of the 
problem based on new precise r-decay data is reasonable. 

The spin-1 parts Hy (q 2 ) and IT^ (q 2 ) are analytical functions in the complex g 2 -plane 
with a cut along the right semiaxes starting from the threshold of the lowest hadronic state: 
4m 2 for Uy^ and 9m 2 for IT^ . The latter has a kinematical pole at q 2 = 0. This is a 
specific feature of QCD, which follows from the chiral symmetry in the limit of massless 
u, <i-quarks and its spontaneous violation. Indeed, in this limit the axial current is conserved 
and there exists a massless Goldstone boson, namely the pion. Its contribution to the axial 
polarization operator is given by: 

n». = fl(^- qj f) , (3) 

where is the pion decay constant, f n = 130.7 MeV [f|. When the quark masses are taken 
into account, then in the first order of quark masses or, what is equivalent, in m 2 , eq. (f|) 
gets modified: 

KM,- ■ (4) 

It can be decomposed in the tensor structures of (jl]): 




According to this equation the residue at the kinematical pole is equal to —f 2 - The accuracy 
of this statement is of order of the chiral symmetry violation in QCD, ~ m 2 /m 2 , where m p is 
characteristic hadronic scale (say, p-meson mass) (e.g. a subtruction term ~ g^f^ml/m 2 
can be added to IT^,). 

The difference HV — II^ is of particular interest, since in QCD it does not have any 
perturbative contribution in the limit of massless quarks. We use the analytical properties 
of Ily^s) and H^\s) in the complex s-plane in order to construct the sum rules for Uy — IL4 
valid at large \s\. At large \s\ the operator product expansion (OPE) takes place 

n« (s) - n« W = £ <gj ( 1 + c D ^) = Y, T^m • < 6 > 

D>4 ^ ' D>4 ^ ' 

where 0]^~ A are the vacuum averages of local operators, constructed from quark and gluon 
field. In what follows the operators Od without index V — A include the radiative corrections 
Od = 0^ A (1 + cdolsIk). Higher order perturbative corrections to 0^~ A , as well as the 
terms ~ m 2 ld are neglected. One may expect, that OPE is valid in the whole complex 
s-plane, except for the domain of small \s\ and near positive real semiaxes (see Fig. |T|). 

The measured difference of the spectral functions v 1 (s) — a\ (s) is shown in Fig. |^. In this 
paper we use the ALEPH data, since the files with invariant mass spectra and correspondent 
covariance matrices are publicly available. 
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Figure 1: Domain of OPE validity 



The operators 0^ A have been computed up to dimension D = 6 in ]7||. The earlier 



calculations of have been done in P|— 1|10||, but there are some discrepancies in the results. 



We have recalculated 0$' A (see Appendix). In the calculation of Oq~ a and 0%~" A the 
factorization hypothesis, i.e. the saturation by intermediate vacuum state, is assumed. As 
shown in Appendix, there is an ambiguity in the factorization of D = 8 operators among 
the terms of order N~ 2 ; they are neglected here. The results are: 

0\~ A = 2(m u + m d ) <qq>= (7) 
0\~ A = 27ra s <(^A a rf)((i7MA a w)-(^757MA a rf)(^757^ a M )> 

= g" i <OT> (8) 

Ol~ A = 8na s m 2 <qq> 2 . (9) 

The definition of itiq is given in Appendix, we assume the isotopic symmetry among the 
quark condensates: < uu >=< dd >=< qq >. 

Let us discuss what is known about the vacuum averages 0^ ) ~ A . The numerical value 
0\~ A = —f^m 2 = —3.4 x 1CT 4 GeV 4 is very small and in almost all cases can be ignored. 
The quark condensate < qq > can be found from Gell-Mann-Oakes-Renner low energy 
theorem 

<qq>= -r J \" (10) 
2 m u + m d 



At standard values (see e.g. |I2|) m u = 4.2 MeV, nid = 7.5 MeV we have 

< qq >= - 1.4 x 1(T 2 GeV 3 (11) 
The value of < qq > depends on the normalization point /i 2 and it is unclear to which normal- 



ization point it refers. In recent analysis of QCD sum rules for proton |T^] the same numerical 
value as ( ]TT| ) was found at the point fx 2 = 1 GeV 2 . Using this value and a s (l GeV 2 ) = 0.5, 
which follows from a(m 2 z ) = 0.119 by using three loop QCD renormalization group evolution, 
we get for renorminvariant quantity: 

a s <q q> 2 = 1.0 x 10~ 4 GeV 6 (12) 
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Figure 2: The measured difference Vi(s) — a\(s). Figures from |CT and 



Here, however, we have to be careful. In QCD sum rule analysis no a s corrections were 
accounted. They may result in 20 — 30% uncertainty in < qq >ic e v 2 - Taking (|l|), we get 
for Ol~ A : 



O 



V-A 
6 



- 2.2 x 10~ 3 GeV 



The value of Wq was found in |14j from the analysis of QCD sum rules for baryons: 



0.8 ± 0.2 GeV 2 



The substitution of flTJ) and (0) in © gives: 



Ol~ A = 2 x 10" 3 GeV 8 



(13) 



(14) 



(15) 



Perturbative a s corrections were calculated for the contribution of D = 4 |15j and D = 6 



TBI, O operators. The correction to 0\ is C4 = 4/3; it increases the effective value of the 



operator O4 on 20%. Concerning Oq, two essentially different values have been obtained: 
cq = 247/48 in |TB[ and cq = 89/48 in |T7J Q In [EH] it was argued, that the last one is more 
reliable, since in its calculation the correct treatment of 75 in dimensional regularization 
scheme was done and more plausible vacuum saturation of 4-quark operators was performed. 
If we take c 6 = 89/48, put a s (l GeV 2 ) = 0.5 and neglect the q 2 dependence, then we get for 
the effective operator: 

6 = - 3 x 10" 3 GeV 6 (16) 

In leading order Og weakly depends on the normalization point. So we will consider Og = 
Og~ A as effective D = 8 operator with a s correction included. 

Our goal is to find Oq and Og from r decay data and compare them with ( |T6"D and ([15]). 
Higher order operators with D > 10 also contribute to OPE (|6|). OPE is an asymptotic 

1 There is also the logarithmical correction (a s /47r) In (s/n 2 ) to the operator Oq, which was not included 
in (JsJ) . However this term is small for physically reasonable values of the scale /j, 2 . 
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series. The comparision of numerical values (|T3|) and (|T5|) indicates, that at \s\ = 1 GeV 2 
this series starts to diverge at D = 8 (the same conclusion \Oq\ ~ |Og| in GeV follows also 
from our final result). Therefore in order to get reliable results we have to go to higher 
\s\ or to improve the convergence of the series. In order to estimate the error in the 06,8 
determination we will accept the conservative assumption, that Od measured in (GeV) D 
increase starting from D = 10, for instance |Oio| ~ 2 1 Oe | - 



2 Moments sum rules 

The dispersion relation or Borel transformation requires the knowledge of the spectral func- 
tions for all s. Although the vector function t>i(s) within isotopic symmetry can be found for 
s > m 2 . from e + e~ annihilation, the precision is low, since the experimental analysis involves 
the states with 6 mesons and more. The axial- vector function ai(s) is not known beyond 
this point at all. 



The technique of the spectral moments |T8[ used for the evaluation of hadronic r-decay 



branching ratio does not need this information. The following moments are computed: 
M kl (so) = ± r-(l--) k (-)\vi-a 1 )(s) (17) 



2?r 2 Jo s o V s oJ \ s o 

Cj are binomial factors (we take O2 = / 2 here) . For sq < m 2 the moments can be computed 
from experimental data. In the equation fll8D the integral goes counterclockwise over the 
circle with radius s . 

In principle one can find all operators Op in this way. Nevertheless, for k < 2 the ex- 
perimental error is very high, so the number of independent moments in flT8|) which can be 
computed with desirable accuracy is less, than the number of unknown operators. Conse- 
quently we have to neglect the contribution of higher dimensional operators, introducing 
thereby a theoretical uncertainty. 

In order to find the operators up to D = 8, one should compute four independent moments 
M kl . The experimental error is large for small k and large /. The theoretical uncertainty 
grows with k + I, since unknown operators from 0\q to 02{k+i+i) are involved. Although 
the experimental error could be in acceptable range, the result depends on particular set of 
moments. 

On the other hand, / 2 and O4 are known from other data with high accuracy. One may 
use this information and moments with k = 2 in order to find the operator of dimension 6 
and higher: 



-V0 2n = si? 



n—3 p2 

Y(n - I - 2)M 2l (s ) + (n - 2) ^ + (n - I) 



Oi 



1=0 



n>3. (19) 



Provided that the OPE (0) works, the r.h.s. of this equation should not depend on s - It 
is plotted versus s in the Fig. Oa,b for n = 3 and n = 4 respectively. According to these 
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Figure 3: Right hand sides of the equation (|T9| ) for n = 3 (a) and n = 4 (b). 

figures the operator Oq can be estimated as —(5 ± 3) x 10 _3 GeV 6 , while the operator 0$ 
is even remotely does not look as a constant. The uncertainty in the determination of f v 
strongly affects the result. In the Figs [5]a,b we have plotted the operators for 3 different 
cases: the central value f n = 130.7 MeV and with ±1.5% excess. 

The reason of this failure is the invalidity of the expansion (P) in all complex s-plane. In 
the moments (|18|) the integral over the circle crosses the area where OPE does not work (see 
Fig. |l|) and it is questionable whether this contribution is suppressed enough by the factor 
(1 — s/so) k . As Fig. |3] demonstrates, this is true only for the radius of the circle greater than 
~ 2 GeV 2 . 

In principle eq. ([19]) can be used for n > 5, but the experimental error in this case is so 
high that it does not allow us to extract any reliable information about the values of D > 10 
operators. 



3 Borel sum rules 

The Borel sum rules can be cor 
unsubtructed dispersion relation 



The Borel sum rules can be considered at complex values of s. Represent 11$ — TL^y via 



n« W -n« W = ^r^iW* + fi (20) 

The last term in the r.h.s. is the contribution of the kinematic pole. Let us substitute the 

OPE (|^) in the l.h.s. of (^). Consider s in the complex plane s = sqc 1 ^ (0 = on the 

upper side of the cut) and perform Borel (Laplace) transformation of (|20|) by s . The real 
and imaginary parts give us the following sum rules: 

exp(— cos^cos(— BB^fr-aM— = £ + ^(-)* 2fc (21) 
' ' k=i 
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exp(^cos0j Bin(— sm^j^-aO^)^^ = £(-) fc , M2fc+2 (22) 



fe=i 



The expression in the exponent is negative for 7r/2 < <p < 3ti/2. Since eq. fl2ll) is symmetric 



and eq. (22) is antisymmetric in the lower half plane, it is enough to analyze the region 

7r/2 < < 7T. 

At certain angles the contribution of some operators vanishes. This fact can be used to 
separate the operators from each other. In particular, eq. (|2Xf) with = 37r/4 and eq. (|22|) 
with = 27r/3 do not contain the operator of dimension 8. For <fi = 5n/6 the operator 



0% disappears from the eq. (|2T|) and mainly the operator 0% contributes to the excess over 
f 2 . All these cases are shown in Figs |]|| We also show eq. fl2"2"|) for <fi = 3tt/A, where the 
contributions of the operators 0$ and Os are comparable. Thin areas on the graphs are just 
because the sin or cos for particular <ft and M 2 has zero at s — m 2 . , where the experimental 
error is high. 

Borel sum rules have serious advantage, since the operators of higher dimensions are 
factorially suppressed. This allows one in the sum rules to go from above up to M 2 w 1 GeV 2 
and even lower in some cases. But they have also a disadvantage: at M 2 > 1 GeV 2 the upper 
tail of the integrals in the l.h.s.'s of (|2T|), (|22| ) are not suppressed enough. But, luckely, the 
oscillating factors in the l.h.s.'s of (21]), (|22"D help in some cases as can be seen from Figs ||.^]. 
We exploit this fact. 

Let us look first at eq. (|2l|) at (p = 5n/6. The r.h.s. of (|2l|) is equal to: 

f2 ,V30 i 10 1 _J_Oio . , 

2 M 2 4 M 4 48 M 8 ' { } 

higher orders are discarded. As seen from Fig |]b at M 2 = 0.8 GeV 2 , the deviation from 
f 2 is definitely outside the limit of errors. The second term in ( PB"D is small ~ —3.0 x 
10 -4 GeV 2 . The main contribution comes from the operator Oq, since 0\q contribution is 
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strongly suppressed. If we neglect it, then one gets from Fig |]b: 

6 = -(6.3 ± 1.4) x 1(T 3 GeV 6 (24) 

Possible contribution of O w (at O w ~ 2|0 6 |) and 3% uncertainty in f 2 increases the error 
to 2.5 x 1(T 3 GeV 6 (the errors are added in quadratures) and finally we get from eq. (|2"l|) at 
4> = 57r/6: 

6 = -(6.3 ± 2.5) x 1(T 3 GeV 6 (25) 

This value can be checked by considering the sum rule (22) at <p = 2tt/3 (Fig|5]a). The r.h.s. 
of (£2p reads: 

_ yl 9± _ yl 9*. yl £™_ (2 6) 

2 M 4 4 M 6 48 M 10 1 j 

The most suitable M 2 is in the region of the isthmus in the experimental errors area, M 2 w 
0.85 GeV 2 . Here, according to Fig |^a, the l.h.s. of is (5.3 ± 1.0) x 10~ 3 and we have 
from Q2"5D 

6 = -(6.8 ± 2.1) x 10~ 3 GeV 6 (27) 

in agreement with (|24] ) (the error from O 10 is included). 

Let us try to find the value of the operator 0§. In ([H]) at <fi = 37r/4 the contribution of 
Oq vanishes and in the r.h.s. we get: 

f 2 + ±£± - J-£i. - 19r (ox) 

n y/2 M 2 Q^2 M 6 24 M 8 1 1 

The most appropriate domain of M 2 is the area of small M 2 , where the deviation from f 2 
is remarkable. At the assumption |Oio| ~ 2 1 06 1 the minimal squared error is achieved at 
M 2 = 0.6 GeV 2 

Os = (6 ± 8) x 10" 3 GeV 8 , (29) 
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which gives us only the upper limit of 0§ < 14 x 10~ 3 GeV 8 . Similar upper limit follows 
from consideration of large M 2 rs 1.4 GeV 2 , where the contribution of 0\q operator is small 
and experimental error dominates. 

Consider finally the eq. (^) at <p = 37r/4, where both operators Oq and Os contribute 
(Fig ^|b) . This value of 4> has the advantage, that 0\q operator disappear from the r.h.s. of 
(f22j), which becomes: 

_ ^(h _ lOe _ _J_Os_ 1 Oi2 , , 

y/2 M 4 2 M 6 Q^2 M 8 120^ M 12 1 ' 

The small numerical factor in front of 0\2 operator allows one to go to low values of M 2 , where 
the experimental errors are small. We choose M 2 = 0.65 GeV 2 . Then, even if Oyi ~ 5|Og|, 
its contribution to fl30|) is small. At M 2 = 0.65 GeV 2 the data give the value (8.5±0.6) x 10" 3 
for the expression (0). The substitution of 6 = -(6.8 ± 2.1) x 10~ 3 GeV 6 given by (|27|) 
results to: 

8 = (7 ±7) x 10~ 3 GeV 8 (31) 

(The possible error from Oi 2 contribution is accounted at | O i2 1 ~ 5|0 6 |, all errors are added 
quadratically.) Again, only the upper limit. 

More definite result for 8 can be obtained if we accept more optimistic assumption, that 
the magnitudes of Oio,i2 operators in GeV are of the same order as \Oq\. Then the error of 
Oq in eq. ( p7[ ) is reduced to 1.6 x 10~ 3 and in eq. ( j22"| ) at <f> = 37r/4 one may go down to 



M 2 = 0.4 GeV 2 to minimize the total error. In this case our best values from Borel sum 
rules are: 

6 = -(6.8 ± 1.6) x 10~ 3 GeV 6 , 8 = (7.2 ± 3.4) x 10~ 3 GeV 8 (32) 

These results must be taken with a certain care, sine the errors may be underestimated: 
at such low M 2 there could be some terms, not given by OPE (e.g. of exponential type, 
exp(-/V=s)). 



4 Gaussian sum rules 

In Borel sum rules the spectral functions are integrated with the weight function e~ s l Ml . 
This exponent suppresses the contribution of the points near s = m 2 with low experimental 
accuracy and unknown tail beyond them. However this suppression is not always enough, 
especially when one would like to find the excess due to operators 06,8 over dominating f 2 . 

Gaussian sum rules have an advantage, that the high energy tail in the dispersion integrals 
are suppressed by the factor e~ s l M , stronger than in Borel sum rules even at M 2 not much 
lower m 2 . But they also have a disadvantage, because the factorial suppression of high order 
terms in OPE starts in fact at operators of very high dimension. 

The sum rules of this kind can be constructed with the help of the analysis of the 
correlators on the complex plane. Consider for instance the real part of the polarization 
operator on the imaginary axes: 

RelW, So ) = — ^ __ ads = g_ (33) 
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Since both sides of this equation depend only on Sg, one may apply the Borel operator over 
this variable to get the following gaussian sum rule: 



2tt 2 



e- s2 l M \v x 



ax){s] 



s ds 



■-YO 



Ik 



M 4 f^{k-l)\M± k 



(34) 



Since the operator O4 is negligible, the expansion in the r.h.s. starts from Os- Consequently 
the eq. (|3~4"D can be used to find the operator of dimension 8. The l.h.s. of ( |3~4"D is plotted in 
Fig. |a. 

In order to find the operator 0§ from Gaussian-like sum rule, one has to construct an 
another function of Sq from the correlator IT(s) which consists of the operators of dimension 
Ak + 2. To kill the leading term f 2 , consider the imaginary part of e* v TI(soe l¥ ') at some angle 
ip. Further, to cancel the operators of dimension 4k, one can add this function at the point 
symmetric with respect to imaginary axes. The result is: 



1 

27 



Im 



3^/2 n 



(1) 

V-A 



s e 



+ e 



/2tt(1) , 



-s Q e 



-i<t>/2\ 



smi 



2tt 2 



(vi - ai)(s) 



S 4 + Sq — 2s 2 Sq cos 



s 2 ds = sin (k(f)) 



O 



ik+2 



-,2k+l 



k=l 



Applying the Borel operator by variable Sq, we get: 



'0 



(35) 



2tt 2 



( s2 



COS ( 



sin 



sin (p) (v\ — ai)(s) — %- 
M 4 / V A ; M 2 



sin (k(p) 



O 



■lk+2 



k=l 



k\ M 4fc + 2 



(36) 



The r.h.s. starts from 6 . For the exponent to be decreasing, the angle <fi must be in 
the range n/2 < <p < tt (the range tt < <p < 3n/2 obviously does not contain any new 
information). The l.h.s. of ( p6| ) for = 3tt/A is plotted versus M 2 on Fig. |5]b. 
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The operators O^, k > 1 contribute to the sum rule (|34]). If we neglect contributions of 
high orders starting from O12, then from the Fig ^|a at M 2 = 1.5 GeV 2 for Og we would have 

Os = (10.6 ± 3.6) x 10~ 3 GeV 8 (37) 

However, the result strongly depends on O12 operator. If we use the same estimation as in 
the previous section | O12 1 = 5 1 Oe | , then (|37D may change on ±7.7 x 10~ 3 . Considering this 
amount as possible error in ( |3"TD and adding the errors in quadratures, we get: 

Os = (10.6 ± 8.5) x 10~ 3 GeV 8 (38) 

Going to lower energies is dangerous, because the contribution of 0\ 2 increases drastically. 
At higher M 2 , where the higher order operator can be neglected, the experimental error does 
not allow to get any definite conclusion about the magnitude of 8 . 

Now we turn to eq. fl36|) at = 3tc/A. The most suitable scale is M 2 = 1.5 GeV 2 . The 
next to Oq in (|36|) is the operator 0\q. Its contribution at 1.5 GeV 2 is suppressed not so 
much, by a factor (v^M 4 ) -1 . If we allow, that | O i0 1 could be as large as 2|0 6 | (in GeV) and 
include this uncertainty as an error, then the following esimation goes from Fig §3: 

O e = - (7.2 ± 5.1) x 10~ 3 GeV 6 (39) 

In case of more optimistic assumption used in previous section |Oio,i2 1 ~ 1 06 1 we get 
better results, especially for the operator 8 : 

6 = -(7.7 ± 3.2) x 10~ 3 GeV 6 , 8 = (9.8 ± 2.3) x 10" 3 GeV 8 (40) 

One may stress, however, that the assumption |0 12 | ~ \O e \ is the most dubious. 

The conclusion is the following. The operators, obtained from Gaussian sum rules are 
compatible with Borel ones. However the range of applicability is different. Indeed, the 
Borel exponent effectively suppresses the high energy contribution only for M 2 < lGeV 2 . 
But in the Borel expansion each operator of dimension D has the factor l/(D/2)\, which 
provides much stronger suppression then the Gaussian factor l/(D/4)!. Consequently the 
effective radius of convergence of borel series could be lower. As the Figs f|-[5] show, this is 
indeed true: the coincidence of the right and left hand sides begins with M 2 = 0.6 GeV 2 , 
twice lower the correspondent gaussian value. 

5 Summary 

The recently obtained data by ALEPH and OPAL collaborations on V — A spectral functions 
in r-decay were used for determination of quark and quark-gluon condensates: vacuum 
expectation values of dimension 6 and 8 operators 0\~ A and 0\~ A . The analytical properties 
of polarization operator Uy\q 2 ) — n^(g 2 ) in the complex q ,2 -plane were exploited. Three 
types of sum rules were used: moments sum rules, Borel and Gaussian ones. The results are 
summarized in Tables 1,2. They are in agreement with one another in the limit of errors 
and the best values of 06,8 are: 

Oq = -(6.8 ± 2.1) x 10~ 3 GeV 6 , 8 = (7 ± 4) x 10" 3 GeV 8 (41) 
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source 


assumption 


scale 


central 


exp. 


3/0 /„- 


L^IO 


total 






in vj ( v 


value of 06 


error 


error 


error 


error 


eq. (|21|) at 


O 10 ~ 20 6 


0.80 


—6.3 


1.4 


1.3 


1.6 


2.5 


(j) = OTT/b 




0.62 


—o.l 


0.7 


0.8 


1.3 


1.7 


eq. (p^) at 


Oio ~ 206 


0.85 


-6.8 


1.4 




1.6 


2.1 


<f> = 2vr/3 


Oio ~ 06 


0.85 


-6.8 


1.4 




0.8 


1.6 


eq. (H) at 


Oio ~ 206 


1.53 


-7.2 


2.8 




4.3 


5.1 


= 3vr/4 


Oio ~ 06 


1.47 


-7.7 


2.0 




2.5 


3.2 



Table 1: Values of the operator Oq with possible errors in 10~ 3 GeV 6 , obtained from Borel 
and Gaussian sum rules. In each case the scale M 2 is choosen in such way, that the total 
squared error (the sum of all squared errors), is minimal. In second column the magnitudes 
of operators are given in GeV. 



source 


assumption 


scale M 2 


central 


exp. 


3% f 2 


06 


010,12 


total 






in GeV 2 


value of Os 


error 


error 


error 


error 


error 


eq. (HD at 


Oio ~ 20 6 


0.65 


6.2 


2.8 


1.2 




7.4 


8.0 


<f> = 3vr/4 


Oio ~ 06 


0.59 


5.5 


1.3 


0.9 




4.1 


4.4 


eq. (||) at 


O12 ~ 506 


0.65 


7.0 


1.0 




5.8 


4.0 


7.1 


(p = 3vr/4 


O12 ~ Oq 


0.41 


7.2 


0.1 




2.8 


2.0 


3.4 


eq. (|D 


0\2 ~ 506 


1.49 


10.6 


3.6 






7.7 


8.5 




O12 ~ Oq 


1.29 


9.8 


1.0 






2.0 


2.3 



Table 2: Values of the operator Os with possible errors in 10 3 GeV 8 . 



The errors here are not quite well defined, they are just our estimations based on the data, 
presented in Tables 1,2. Particularly, in case of 0§ operator the errors strongly depend on 
the assumption about the magnitude of Oi 2 . In the most pessimistic case of large 1 12 1 , say 
| O12 1 ~ 5|0 6 | in GeV, we have only the upper limit O s < 14 x 10~ 3 GeV 8 . 

The values ( [41~D are by a factor 1.5 — 2 larger, than the values (|T3D , (|15"D found from low 
energy theorems and QCD sum rules (see Introduction). If this discrepancy is addressed 
to the quark condensate, then, in accord with fllCf ) it means that m u + m d at 1 GeV 2 is by 
20—40% less than the standard value 12 MeV. Up to this may be not so essential discrepancy 
the analysis of r-decay data confirms the general concept of OPE and the magnitudes of 
quark and quark-gluon condensates. 

The authors are thankful to A.Oganesian for his help in the calculation of dimension 8 
operator. The research described in this publication was made possible in part by Award No 
RP2-2247 of U.S. Civilian Research and Development Foundation for Independent States of 
Former Soviet Union (CRDF) and by Russian Found of Basic Research grant 00-02-17808. 
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Appendix: The condensate of dimension 8Q 

It consists of three different contributions: 

8 = of + of + Of , (42) 

where the upper index denotes the number of quarks in vacuum. The purely gluonic conden- 
sate Of and two-quark one Of have been computed in ||. They contain many independent 
operators, which cannot be expressed in terms of condensates of lower dimensions. However 
in the masseless quark limit these operators are the same for vector and axial correlators 
and disappear in the difference 0%~ A . 

We have explicitly computed the four quark condensate for the vector current correlator: 

2 

O g V(4) = 9 -(%{u la \ a D p d){d la \ a Dpu) - llgf abc G c af3 (u lal5 X a d)(d 1/3l5 X b u) 

-14 (uD alf3 j 5 X a ^ Q d)(d lpl5 X a u) - U(dfi a ^ 5 X a ^ a u)(u W X a d) 

+ (u la {G al3 , X a }d)(d ll3l5 X a u) + (d la {G aP , X a }u)(u ll3l5 X a d) ) , (43) 

where G a p = § X a G a a(3 is gluon field, G a/3 = \e a p^G^ v , e i23 = 1, 7s = «7o7i72 7 3, D= l5-*D, 
Dfj, = — l -^-X a A a ^ the derivatives in the brackets (. . .) act only on the objects inside these 
brackets; A a are Gell-Mann matrices, [A a , X b ] = if abc X c . One may check, that the eq. (|43|) 
can be brought to the form obtained in |J, which verifies our results. 

The condensate of axial currents 0% ^ can be easily obtained from ( f4*3"|) with the help of 
the replacement d —>■ ^y 5 d. 

To reduce the number of independent operators in (|43|) , the vacuum insertion can be 
applied to Of . Nevertheless this procedure is not unambiguous. Indeed, let us consider the 
following operator, which appears in the derivation of fl43|): 

g 

4 



< (q la X a q)D (a D p) (q lp X a q) > = § < f abc G c ap {q la X a q){q lp X b q) > 



= - \ < qGq >< qq > , (44) 

where G = ^afiGap. In (|4]) we used the quark equation of motion and commutational 
relation for the covariant derivatives, and then applied the vacuum insertion. On the other 
hand, one may apply the vacuum insertion at first and use equations of motion after then: 

< (qiaX a q)D {a D /3) (qjf 3 X a q) > = - (l - -^j < qD 2 q >< qq > 

= -7T ( 1 ~ T7t I < >< M > > ( 45 ) 



2 V N*, 

where N c is the number of colors. We see, that two different ways of the vacuum insertion 
give the same result only up to the terms of order ~ V-^c- 

2 A.Oganesian participated in the calculation of dimension 8 operator. 
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Consequently, within the framework of the factorization hypothesis we may write the 
four quark operators ( pE3"D in the following form: 

Ol (4) = - Ot {4) = Ana s m 2 < qq > 2 ( 1 + 0(N~ 2 ) ) (46) 

The parameter m is introduced as < qGq >= im^ < qq > according to [T(J. Within this 
accuracy the equation (|46|) coincides with the result of ref. ]10| (according to the isotopic 
symmetry the current correlators Il^ erc ' ) = 2Il^ f ' 0- 1 ). 
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